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Abstract: We consider Langevin equation involving fractional Brownian 
motion with Hurst index H 6 (0,^). Its solution is the fractional Ornstein- 
Uhlenbeck process and with unknown drift parameter 6. We construct 
the estimator that is similar in form to maximum likelihood estimator for 
Langevin equation with standard Brownian motion. Observations are dis¬ 
crete in time. It is assumed that the interval between observations is n —1 , 
i.e. tends to zero (high frequency data) and the number of observations 
increases to infinity as n m with m > 1. It is proved that for positive 0 
the estimator is strongly consistent for any m > 1 and for negative 0 it is 
consistent when m > . 
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1. Introduction 


Let (12, S’, P) be a complete probability space. We consider fractional Brownian 
motion B H = {Bp , t > 0} on this probability space, that is, the centered 
Gaussian process with the covariance function 

«(t, »)=i («“+<“-1*-*n- 

We restrict ourselves to the case H £ (0, and consider the continuous (and 
even Holder up to order H) modification that exists due to the Kolmogorov 
theorem. Introduce the Langevin equation, 

X t =x 0 + e [ X.ds + B?, t > 0, He( 0,±). (1) 

Jo 


According to Proposition A.l from [7], this equation has the unique solution 
that is named fractional Ornstein-Uhlenbeck process and can be presented as 

X t = x 0 e et + de et [ e~ 9s BP ds + B?, t > 0. (2) 

Jo 

The goal of the paper is to construct consistent (strongly consistent) estima¬ 
tor of the unknown drift parameter 6 by discrete observations of the process 
A'. 

The problem of the estimation of the drift parameter 6 in the linear equation 
containing fBm and in the equation (1) when the Hurst index H > ^ was 
investigated in many works. For linear models, mention only papers [2] and 
[13]. Drift parameter estimators for fractional Ornstein Uhlenbeck process with 
continuous time when the whole trajectory of X is observed, were studied in 
[1, 12, 15]. Kleptsyna and Le Breton [15] constructed the maximum likelihood 
estimator and proved its strong consistency for any 0 6l. They also investigated 
the asymptotic behaviour of the bias and mean square error of this estimator. 
The sequential maximum likelihood estimation was considered in [20]. Hu and 
Nualart [12] proved that in the ergodic case (9 < 0) the least square estimator 


0t 


So XtdX t 

So X tdt ’ 


( 3 ) 


is strongly consistent for all H > | and asymptotically normal for H £ [i, |). 
They also obtained the strong consistency and asymptotic normality of the 
estimator 

^=(-HmfJS x ' dt ) " (4) 

In [1] the corresponding non-ergodic case 9 > 0 was considered and the strong 
consistency of the least square estimator (3) was proved for H > J. It was 
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obtained also that e 9t yO t — Oj converges in law to 20C(1 ) as t —> oo, where C(l) 
is the standard Cauchy distribution. Minimum contrast estimators in continuous 
and discrete case were studied in [4] . The distributional properties of maximum 
likelihood, minimum contrast and least square estimators were explored in [21]. 
For the two-parameter generalization see [8] . 

In [6, 10, 11] the discretized version of (3) is considered, namely 


Or,. = 




( 5 ) 


where the process X was observed in the points f•; = iA n , i = 0,... ,n, such 
that A„ — > 0 and nA n —> oo as n —> oo. In [6] the ergodic case 0 < 0 was 
studied, the strong consistency of this estimator was proved for H > | and the 
almost sure central limit theorem was obtained for H e (i, |). The non-ergodic 
case 6 > 0 was considered in Es-Sebaiy and Ndiaye [11]. They proved the strong 
consistency of the estimator (5) for H e (|,l) assuming that A„ —> 0 and 
—> oo as n —>• oo for some a > 0. The same result was obtained for the 
estimator 


0 n = 


X l 


2a„e”=i x i _: 


In [14, 23] the following discretized version of the estimator (4) was considered 


0 n = — 


nHT(2H) 


k-1 / 


where 0 < 0 and the process X was observed in the points A,2A, ...,nA 
for some fixed A > 0. Hu and Song [14] proved the strong consistency of the 
estimator for H > \ and the asymptotic normality for \ < H < |. 

In [5, 24] more general situation is studied, where the equation has the 
following form dX t = 0X t dt + adBf 1 , t > 0, and d = ( 0,a,H ) is the un¬ 
known parameter, 0 < 0. Consistent and asymptotically Gaussian estimators of 
the parameter 0 are proposed by the discrete observations of the sample path 
(XkA n ,k = 0,..., n) for H e (^, |), where nA^ —> oo, p > 1, and A„ —>• 0 as 
n —> oo. In [24] the strongly consistent estimator is constructed for the scheme, 
when H > i, the time interval [0,T] is fixed and the process is observed at 
points h n , 2 h n ,..., nh n , where h n = E 

In [9,17] the so called sub-fractional Ornstein Uhlenbeck process was studied, 
where the process B^ in (1) was replaced with a sub-fractional Brownian mo¬ 
tion. In [9] the maximum likelihood estimator for such process was constructed, 
in [17] the estimator (3) was investigated in the case 0 > 0. The maximum likeli¬ 
hood drift parameter estimators for fractional Ornstein Uhlenbeck process and 
even more general processes involving fBm with Hurst index from the whole 
interval (0,1) were constructed and studied in [22]. These estimators involve 
singular kernels therefore are more complicated to study and simulate. To the 
best of our knowledge, it is the only paper when discretized estimates of the 
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drift parameter are constructed in the case H < \. However, the observations 
of the real financial markets demonstrate that the Hurst index often falls below 
the level of taking values around 0.45-0.49 ([3]). In order to consider the case 
of H < i and to overcome the technical difficulties connected with singular 
kernels, we construct comparatively simple estimator that is similar in form to 
the maximum likelihood estimator for Langevin equation with standard Brow¬ 
nian motion. Observations are assumed to be discrete in time and we assume 
that the interval between observations is n _1 , i.e. tends to zero, so we consider 
high frequency data. At the same time, the number of observations increases 
to infinity with the speed n m with m > 1. Let n > 1, tk, n = 0 < k < n m , 

where m £ N be some fixed number. Suppose that we observe X at the points 
{tk, n ,n > 1,0 < k < n m }. Consider the estimator 


6 n {m) 


_„ m 1 

ELo Afc [7I AA'fc |Tt 


_1 v^n m -l y2 
n 2-^k—0 k,n 


( 6 ) 


where X — Ab fc n , XXf^ n — X k.n ■ 

By (1), estimator 9 n (m) from ( 6 ) can be represented in the following form, 
which is more convenient for evaluation: 


0 n (rn) 


°T,k=0 ^k,njk n (A s Xk,n) ds + Efc=0 -^-k,nXB^ n 

1 V2 

n 2-^/k —0 k,n 


( 7 ) 


It is proved that for positive 6 the estimator is strongly consistent for any m > 1 
and for negative 9 it is consistent for m > 277 . 

Our paper is organized as follows. In Section 2 we consider auxiliary result, 
namely, bounds with probability 1 for the values and increments of fractional 
Brownian motion and fractional Ornstein Uhlenbeck process. The bounds are 
factorized to the increasing non-random function and random variable not de¬ 
pending on time. In Section 3 we get the bounds for the numerator of the estima¬ 
tor, while in Section 4 we relate discretized integral sum in the denominator of 
the estimator to the corresponding integral f Q X^ ds. This relation is convenient 
for some values of parameters because it is easier to apply L’Hopital’s rule to 
the integral f Xf ds than Stolz-Cesaro theorem to the sum Efc ^ 1 n 
terms depending on n. Section 5 contains two main theorems, Theorem 5.1 and 
Theorem 5.4 on strong consistency for 9 > 0 and consistency for 9 < 0. Section 
6 contains some auxiliary results and Section 7 is devoted to numerics. 


2. Bounds for the values and increments of fractional Brownian 
motion and fractional Ornstein Uhlenbeck process 

In what follows we shall use next auxiliary estimates for the rate of asymptotic 
growth with probability 1 of the fractional Brownian motion and its increments. 
Throughout the paper considering functions of the form f p log t, p > 0 we 
suppose that 0 ■ 00 = 0 . 
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Proposition 2.1. (i) For any p > 1 and any H £ (0,1) there exists nonneg¬ 

ative random variable £(p,H) such that 


sup |5f | < ((i*|logi|P) Vl)f(p,JT), (8) 

0 <s<t 


and there exists such number c^(p, H) > 0 that for any 0 < y < c^(p,H), 
Eexp{y£ 2 (p, H)} < oo. 

(ii) For any q> \ and any Ff £ (0,1) there exists nonnegative random variable 
rj(q, Ff) such that for any 0 < t\ < £2 < 00 

\ B t*~ B t x \ < (|log(t 2 — ti)| 1/2 + l) (\og(t 2 + 2)) q r 1 (q 1 H), (9) 


and there exists such number c v (q 1 Ff) > 0 that for any 0 < y < c v (q, Ff), 
Eexp {ytf(q,H)} < 00 . 

Proof. The 1st statement was proved in the paper [16]. The 2nd statement fol¬ 
lows immediately from the next relation that can be proved similarly to Theorem 
1 from [18], where even more complicated functional than the increment of frac¬ 
tional Brownian motion, more precisely, fractional derivative, was considered. 
So, we have from Theorem 1, [18], that for any q > \ and any Ff £ (0,1) the 
random variable 


v(q,H) 


sup 



tdH tdH I 

\ B C ~ B t 2 \ 


(t 2 -t\) H \ 

(|l0g(*2 

T —1 

+ 

1 

1 (log(t 2 + 2))9 


is finite almost surely, whence (ii) follows. □ 

Now our goal is to estimate the numerator in (7) and compare it to the de¬ 
nominator. At first, we write the bounds for the values of A' and its increments. 


Lemma 2.2. We have the following bounds for fractional Ornstein-Uhlenbeck 
process X in terms of supporting fractional Brownian motion: 


i) Let 9 > 0. Then for any t > 0 

sup \X S \ <\xo\e 6t + 8e et f e~ 6s 
0<s<£ Jo 


r „ sup | B (f | ds + sup | B^ | 

0 0 <u<s 0<s<£ 


and for any s £ [£, *±i) 

sup \X U -X k , n \< f (e eu (\x 0 \+9 [ e- 1 

J -- V V Jo 


sup \B^\ dv 

)<z<v ) 


- I A 

0<z< 

H 

'u -^k^n 


sup | B’J j J du + sup | Bh - B^ 

0 <z<u J i<u< s 

Let 9 < 0. Then for any t > 0 


( 10 ) 


( 11 ) 


( 12 ) 
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and for any s G [|, 


sup \X U 

£<u<s 


Xk,n | < 


I 61 ! kol 

n 


m 

n 


sup 

0<u<s 



+ sup \B*-Bg n 

±<u<s 


(13) 


Proof. (i ) Bound (10) follows immediately from (2), and bound (11) follows 
immediately from (10) and (1). 

(ii) Bound (12) follows from (2): 

\X t \ <\xo\e et + \6\e et sup \B^\■ f e~ Bs ds + \B^\ < |xo| + 2 sup \B^\. 

0<s<£ JO 0<s<£ 

To establish bound (13), we substitute (12) into the following inequality that 
can be easily obtained from (1): for s > ^ 

\Xs - x k , n \ < \e\ £ \X U \du + \Bf - B* n I. 

□ 

Remark 1. Plugging p = 2 and q = 1 m<o </ie formulae (8) -(9), we get the 
following bounds: 

sup |Rf| < (t ff log 2 t + l)£(2,R), (14) 

0<s<£ 


and for s €[$,*£] 

< (s - £)" (|!og (s - |) |' 1/2 + l) log(s + 2 ) 77 ( 1 , H) 


Rf -B? 


< 


((« - X) H |log (s - £) 1 1/2 + (a - £)") log(s + 2 ) 77 ( 1 , H). 


(15) 


Function f{x) = ^llogicl 2 is bounded on the interval (0,1] for any r > 0. 
Therefore 


(s-v H H(s-±)r<c(s^) 


k \|V2 




for any 0 < r < H. Furthermore, for s £ [—, we have that (s — < 

(s — —) . Therefore, we get from (15) that for any 0 < r < H and for 
s G I"-, 

L n ’ n J 


B? -B? 


<C(s- X) H r log(n m_1 + 2 ) 77 ( 1 , H). 

It follows immediately from (14) that for 6 > 0 

[ e~ 9s sup \B*\ds<£{2,H) [ e~ Bs (s H log 2 s + l) ds < Cf(2, H), (17) 
JO 0 <u<s Jo 


(16) 
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and therefore both integrals jf > e~ 0s Bffds and / 0 °° e~ 9s sup 0<u<iS \Bff\ ds exist 
with probability 1 and admit the same upper bound Cf{2,H). Combining (10)- 
(13), (14) and (16), we get that for 9 > 0 

sup \X U \ < \x 0 \e 9s + C9e 0s f(2, H) + (s H log 2 s+ l) £(2, H), 

0 <u<s 

and for se 

sup \X U -X k ,n\ <0 

£<u<s 

+ {u H log 2 it + 1)£(2, H)) du + (n~ H+r \ogn) rj(l,H), 
while for 6 < 0 

sup \X U \ < |x 0 | + 2 (s H log 2 s + 1) £(2 ,H), 

0<u<s 


j k (e du (\xo\+Cei;(2,H)) 


and for s€[^^ i] 


sup \X U - X kjn \ < ^ ^ sup \B%\ 

Tt Tl 0<u<s 


“ <ZL<S 


H nH | ^ l»IW + log 2 s + 1)^(2,^) 

n n v 7 


sup \B U -B kn \ < 

±<u<s 


+ (n H+r logn) 77 ( 1 , H). 


To simplify the notations, we denote by C any constant whose value is not 
important for our bounds. Furthermore, we denote by 3 the class of nonnegative 
random variables with the property: there exists C > 0 not depending on n 
such that Eexpja^ 2 } < 00 for any 0 < x < C. For example, £(2 ,H) + C and 
77 ( 1 , Ft) + C, Cf(2, H ) and Crj(\, H) for any constant C belong to 3- Also, note 
that for fixed m > 1 and n > 3 we have the upper bound log ( n m ~ 1 + 3) < 
Clogn. Moreover, for any a > 0 there exists such n(a) that for n > n(a) we 
have logn < n a . Taking this into account and using the simplified notations, 
we get the bounds with the same £ £ 3- for 9 > 0 we have for any fixed a > 0, 
starting with n > n(a): 

sup \X U \ < (e 9s + s H log 2 s) £ ( 18 ) 

0 <u<s K ' 


and for s £ [£, *± 1 ] 

sup \X U - Xk, n \ < (£ e 9s + ±s H log 2 s + n~ H+a ) £, (19) 

£<u<s 


while for 9 < 0 


and for s £ 


k fc+1 1 
n ’ n J 


sup \X U \ < (1 + s H log 2 s) C 

0<u<s 


sup \X U -X kt n\ < + ^S H log 2 s + n H+a )C 

— <u< s 


( 20 ) 


( 21 ) 
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3. Bounds for the numerator of the estimator 

Now we are in position to bound both terms in the numerator of the right-hand 
side of (7). At first, give the bound with probability 1 for the 1st term in the 
numerator of (7). All inequalities claimed in Lemma 3.1 hold for any a > 0 
starting with some nonrandom number n(a). 

Lemma 3.1. (i) Let 9 > 0. Then for any m > 1 there exists such f € 3 that 



(ii) Let 9 < 0. Then we have two cases. 

(a) Let 1 < m < jj. Then there exists such ( € 3 that 



•2 2H— 1+ck 


(b) Let m> -jj. Then there exists such ( € 3 that 



Proof, (i) It follows immediately from (18) that 



( 22 ) 


n 


Now we take into account (22), substitute instead of s into (19) and apply 
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Lemma 6.1 to get for any a > 0 the following relations: 


n — 1 


^ ^ Xk,n I (X s Xk,n ) ds 


k =0 


n — 1 


c 2 - y 

k—0 

x (^ e ^ + ^(^) Hlog2 (^ i ) +n ” ff+Qlogn ) 

/ , n m -l n m — l 

=c 4 e + 4 e e '“ (^) i°s 2 (^) 


fe =0 
n m — 1 


k=0 


+3 E ( i ?)“lo8 < ( i ?)+’-- I -'' + “log' 


fc =0 
/n rn — 1 


n m — 1 


E«’“+E m'V 24 ? 


k—0 


k—0 


<C 


0 2 


_07i n 


-(■ 


n H(m—l)+m—2+a _|_ n ~H+a' 


+n 


m—l — H+(m—l)H+a 


(23) 

Evidently, the term ^e 26n dominates and other terms are negligible, whence 
the proof of (i) follows. 

(ii) According to (20), now 

I AVI < sup \X U \ < (l + (^iflog 2 *±i) C- 

0<u<ii±i V ' 

Substituting instead of s into (21), we get the following relations: 


71 — 1 


^ ^ Xk,n I {Xs Xk ?n ) ds 


k—0 


n — 1 


< 


10 —j- 

c 2 l E + 




fc=0 


x(t + K i t i )''lo8 2 ( i t i )+»-" + “) 


n — 1 


= C 2 W A (w)"l08 2 (“) + ^ E 


1 "v-^ 1 2H (^4) 

1 ^ x fc+1 V T 4 fc+1 

\ n ) ® 77 , 


k=0 


k =0 


71—1 




E (^)"i»s 2 “ 


k—0 


Substituting the bounds from Lemma 6.1 into the right-hand side of (24), we 
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obtain 

< C 2 ( 'n m ~ 2 + n (m - 1)H - 2+m log 2 n 

+n 2H(m-l)-2+m log 4 ^ + n m-l-H+ a + n (m-2)H-l+m+ a log 2 \ _ 

We take into account that logn = o(n“) as n —>■ oo, for any a > 0. So, it is 
necessary to compare exponents m — 2, (m — 1)H — 2 + m, 2 H (m — 1) — 2 + m, 
m — 1 — H and (m — 2 )H — 1 + to. We get that exponent 2 H(m — 1) — 2 + to 
is the biggest under the condition m > while exponent (to — 2)H — 1 + m is 
the biggest under the condition m < jj whence the proof of (ii) follows. □ 

Now we establish moment bounds for the 2nd term in the numerator of the 
right-hand side of (7). In order to do this, we apply well-known Isserlis’ formula 
for calculation of higher moments of Gaussian distribution: let {xi, X' 2 > X3; X4} 
be a Gaussian vector, then 

E(XlX2X3X4) = E(X1X2)E(X3X4) + E(X1X3)E(X2X4) + E(X1X4)E(X2X3). 

Therefore, we can calculate mathematical expectations EB?? B(? ABj* ABj 1 for 
k ^ j as 

ESf B* ABf ABf = EBfABfEBfABf + EBf ABf EBf ABf 

+ EBf BfEABf ABf < EBf ABf EBf ABf 

+ EBfABfEBfABf (25) 

because for H G (0, |) the increments of fBm B ff are negatively correlated and 
so EABf ABf < 0. Similarly, 

EBf Bf (ABf„) 2 = 2EBf ABf n EBf ABf„ + n" 2 *EBf Bf. (26) 

Lemma 3.2. (?) Let 0 > 0. TTien /or any to > 1 we have the following moment 
bound 

( n m -1 

E 

k =0 

(ii) Let 6 < 0. TTien /or any to > 1 we /lane t/ie following moment bound 

( n rn — 1 

E 

fc=0 

Proof, (i) It follows from (2) that 


* ABj? 


< Cn 2m ~ AH . 


-abL 


< Cn 


2—AH 26n n 


a — 1 

E 

k—0 


X k . 


(X$ 


n m —1 


n m —1 


E Xk,nAB£ n =x 0 E 


fc=0 


k—0 


n — 1 

' E < 

fc=0 


n —1 


_—0 s tjH 


B j* ds • A B^ n + ^ B^ n AB^ n —: J* + 


fc=0 
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Remind that EABf n AB^ n < 0 for k ^ j. Therefore 


/n m — 1 


0 - E E e ^ AB k, 


J k,n L ^ 1J j,n 

\ 2 
. H \ 


k =0 


< n 


-2H 


E 

/c=0 


e 2 ^ < n 1 - 2 ^ 


e 2es ds = Cn 1-211 e 29n ” 


So, E (4) 2 < dn 1 - 2H e 29n ” x , n > 1. 

Consider 7®. It is well known (see, e. g., the relation (1.8) from [19]) that for 
H e (0, \) 


n— 1 


n 2H-l 


J2B£ n AB£ L 


2 1 


fc =0 


Therefore, there exists C > 0 such that 


n—1 


E <<?> »>1- 


k—0 


Now we can use the self-similarity property of B H , namely, 

(R o V>0)^a ff (Bf,t>0) 

and get 


e(j„ 3 ) 2 = e £ B" n AB" 
V fc=0 

/ n m —1 

= n 4(1_m)H E I y] B 1 ^ ( 


r>H r>H 
& fc+i — -D fc 


fc =0 


n m —1 


< n 


At last, 


/ IL — _L 

, n 2 m (l-2H) E I n y gH (gH _ gH \ 

\ ' ^ tT 7717 V n m n Trr / 


k =0 


^ -AmH-\-2m—AmH _ ^j^AH-\-2m—SmH 


n m — 1 

2 „o ^ fi A +e i_ /" r e Ss-8u EB H B H AB H AB H duds 


0 < E (7 2 ) 2 = 0 2 ^ e 

j,k=0 
n m —1 „A k 

20# / n I 


0 Jo 


= ® 2 E 


/c—0 


JO JO 


< »■ *"Kn;'n" (AB« n fduds 


+« 3 E 




e"® 8 " 0 " EB?BZAB? n AB? du ds =: Jf + J 2 ". 




/o Jo 


( 27 ) 
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We get from (25) and (26) that 
,2 1 


E B?B? 


»((*?)“-1--©“+ 

< Cn~ 4H + Cn~ 2H ( s 2H + 


fe+1 


,2 H 


l h |2 H 
\S ~ - 


, u 2H_l s _ u 


/ u \ 2H | l , 1 2H\ 

- (-) +\u--\ 

\n J I n I / 

- n~ 

2H ), 




(28) 


and 


EBfABf< 


fc+i 


+ u — 


((“)“-©“-I- 

- | s - + | s _ L| 2H ) < cfc- 4 * 


i2.ff> 


X (<^) 2 " - (I) 2 " -I 

Substituting the above bounds into (27), we get that 


(29) 


J" < Cn}~ 2H e 26nm 1 and J 2 " < Cn 2 - 4i? e 2e ” m \ 

Evidently, for m > 1 the biggest contribution is from the bound Cn 2 ~ 4H e 20 "™ 
whence the proof follows. 

(ri) Let 9 < 0. In this case 




/n m -l > 

2 n m -l 

0 < E £ e^AI? fc ff n 

S E 

V fe=0 ) 



< 


n m n~ 2H = n m ~ 2H . 


k =0 


So, E(J 4 ) 2 <Cn m ~ 2H , n> 1. 

The term Z? is estimated as before, and 


n — 1 


0 < E (I 2 ) =9 2 Y. e 

j,k =o 

n m — 1 


0£+0- 


_ — 6s—6u Tjy r>H tdH a r?// 


/0 ^0 


EB«B^AB^ n AB" n duds 


=« 2 £ 


„20£ 


/c—0 


'0 ^0 


■'»> dud. 


+« 2 I> 




~-e.-0u EB H B H AB H AB H duds = . J™ + jn 


'0 ^0 


''s L ^ ±J k,n L ^ 1J j,n 


( 30 ) 
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Substituting bounds (28) and (29) into (30), we get that 


to 


e~ es - 0u dsdu 


n m —1 

k k 

r — r — 

Y e 2S " 

/ / e- 9 s - 9 u (s 2H + u 2 H -\s-u 

k=0 

Jo Jo 


n m —1 k 

< Cn 2rn ~ 

±H + Cn -2H £ e 9 ™ r e~ 9s s 2H ds 


fc=0 "'° 


|2 H 


) ds du 


<Cn 2m ~ iH + Cn 1 ~ 2H ( e 9s [ e~ eu u 2H duds 
Jo Jo 

rn ™- 1 M 

<Cn 2m ~ AH + Cn 1 ~ 2H e enm 1 / 


/ e 9u u 2H du + Cn 1 2if / u 

Jo Jo 

/»n m —1 

< Cn 2m " 4fl + Cn}~ 2H / 

Jo 


2H du 


u 2H du 


< Cn 2m ~ 4H + Cn 2Hrn+m ~ 4H < Cn 2m ~ 4H , and J 2 n < CVi 2 " 1-474 . 

Comparing exponents 2m — AH, AH + 2m — 8 Hm and m — 2 H, we get that for 
m > 1 2m — AH is the biggest one, whence the proof follows. □ 

Corollary 3.3. (i) Let 9 > 0. Then 


n m — 1 


E n 4H - 2 e- 29nT ‘ 


1 J2 X *,nAB* n < C. 


k—0 

If we denote f n = n 2H ~ 1 e~ 9nm 1 J2k= o" 1 X k,nABjf n then sup,^ E£ 2 < oo. It 
means that for any m > 1 the numerator of (7) can be bounded by the sum 


c 


■2 -l„20n" 


1-2if 0n" 


f,ni 


where sup ra>1 E£ 2 < oo. 

(ii) Let 9 < 0. Then we have two cases. 

(а) Let 1 < m < jj. Then for any a > 0 the numerator of (7) can be bounded 
by the sum 

£2 n (m-2)ff+m-l+ a + n m-2H ^ 

where sup n>1 E£ 2 < oo. 

(б) Let m > -g. Then for any a > 0 the numerator of (7) can be bounded by 
the sum 

Q2 n {2H+l)m-2H-2+a _|_ n m-2H ^ 

where sup n>1 E£ 2 < oo. 
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4. How to deal with the denominator 


Now our goal is to present the denominator of (7) in more convenient form. At 
first we compare the sum X% n to th e corresponding integral 

f 0 Xf ds. The reason to replace the sum with the corresponding integral 
is that for some values of H and to we can prove the consistency with the help 
of some kind of L’Hopital’s rule, however, the application of L’Hopital’s rule or 
Stolz-Cesaro theorem to the sum X^o" 1 ^kn * s problematic because not only 
the upper bound but the terms in the sum depend on n. 

Lemma 4.1. (i) Let 9 > 0. Then there exists such (3 G 3 that 



X: 


1 


n rn — l 

£ x l 

k =0 


< 


Cl g 2flln m_1 
n 


( ii) Let 9 < 0. Then we have two cases. 

(a) Let 1 < to < jj. Then there exists such Ci G 3 that for any /3 > 0 we have 
the following bound 



1 


l -i 


X 1 ds ~- E X ln 

k =0 


< £-2^mH+m—2H— 1+/3 


(b) Let to > jj. Then there exists such Ci G 3 that for any (3 > 0 we have the 
following bound 



1 n — 1 

xUs-- e x l 

k =0 


n 


^ /-2 2mH-\-m—2H-2-\-(3 


Proof. Evidently, the difference between the integral and the corresponding in¬ 
tegral sum can be bounded as 



X: ds — 


71—1 

E 

k =0 


X, 


k,n 


< / Iv5 n (s)|ds, 

Jo 


where the integrand has the form 

<Pn(3) = {Xl - Xln) 0 < k < U m - 1. 

Furthermore, the integrand can be bounded as 


\(p n {s)\ E \Xs -^Oc,n| (\Xs\ + |^fc,n|) [ k y fc+1 j 

E 2 \X S -Efc,n| sup |-En| Hiy 

0<n<S Ln’ n ; 
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(i) Let 9 > 0. Then from (18), (19) and similarly to (23), 

\X S - X fc ,„| sup |A' U | 1 u k+i\ < (n~ l e 29s + 2n~ 1 e 9s s 11 log 2 s 

L"’ n / V 


+e 9s n H+a + n 1 s 2H log 4 s + n H+r s H log 2 s^Ci • 


(31) 


Integrating over [0, n m_1 ], we see that the integral of the first term in the right- 
hand side of (31) dominates, whence the proof follows. 

(m) Let 9 < 0. Then according to (20)-(21), 


\x.-x k , n \ sup \x u \ i j£ [i mi) < (^ + ^ ) H log 2 (ffr ) 

0<M<S L " " ' ' 

+n~ H+a ) (l+(^i)"log 2 (^i))c i 2 , 


therefore 


k =o 


do \ n I -n 

1 n m — 1 

E (^ i ) 2ffl °g 

n m —1 


H loo ' 2 fc + 1 


l 2ff loo - 4 h±l q. 


fe =0 


fc =0 


(32) 


To get rid of logarithms, we apply Lemma 6.1 to (32) and obtain that for any 

/? > 0 


1 

K(S)| ds < (n m ~ 2 + n mH+m-H- 2+/3 + n m-H-1+0 
_|_ n 2Hm+m-2H-2+/3 _|_ ^ 

Comparing the exponents to — 2, m.H + m — H — 2, m — H — 1, 2Hm-\-m — 2H—2 
and mH+m—2H—l 1 we deduce that for 1 < to < jj the biggest exponent equals 
mH+m—2H — 1, and for to > jj the biggest exponent equals 2Hm+m—2H — 2, 
whence the proof follows. □ 

Corollary 4.2. (i) Let 9 > 0. Then there exists such Ci £ 3 that 


mH-\-m—2H—l-\-/3' s j ^-2 



1 

n 


n m — 1 


E = 



X 2 ds + tin. 


where 


|i?n| < -e 29r 
n 


(ii) Let 0 < 0. Then we have two cases. 
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(a) Let 1 < m < jj. Then there exists such Ci £ 3 that for any /3 > 0 


1 

n 


n m — 1 


E X ln = 



X 2 ds + d n (l3), 


where 


\d n (/3)\ < Cin mH+m ~ 2H - 1+ P. 


( b ) Let m > jj. Then for any f3 > 0 representation (33) holds with 


(33) 

(34) 


\d n {/3 )| < Ci n 2mH+m - 2H - 2+ P. 


5. Main consistency results 


Consider separately cases 9 > 0 and 9 < 0. The case 9 > 0 is more simple and 
additionally the estimator is strongly consistent. The case 9 < 0 needs some 
additional calculations and we prove only that the estimator is consistent. 

Theorem 5.1. Let 9 > 0. Then for any m > 1 estimator 9 n (m) is strongly 
consistent. 

Proof. According to Corollaries 3.3 and 4.2, it is sufficient to prove that 


■2~. —1 „20n TT 


i’n ■ = 


C 2 n 1 e 


,l-2£r„0n” 


& 


Jo" " X 2 ds + d r , 


a.s. as n —> oo, where sup n>1 E£ 2 < oo and |$ n | < A- e 2e " m 1 . Rewrite if n as 

C 2 n -l +n l-2H e - S n™-^ n 


ifn - — 


jn ^2 dg + g-20n— 


Evidently, C^n 1 —> 0 a.s. as n —> oo. Furthermore, for any x > 0 




< 




< 


C 


— x 2 n iH-2 e 28n m - 1 — z .2 n 4ff-2g2en m - 1 ’ 


and the series n 4 H- 2 ^ 2 Sn^-i converges. It means by Borel-Cantelli lemma 

that n 1 ~ 2H e~ enm —»■ 0 a.s. as n —> oo. Evidently, e -20 "™ —>■ 0 

a.s. as n —> oo. At last, according to (2) and the L’Hopital’s rule, 


lim 

T —YCG 


Jo X 2 ds 
e 28T 


X 2 


t^o 2 9e 2BT 


(29) 


-1 


XQ + 


_— 6s tdH 


Bf ds 


and the limit random variable is positive a.s. as the square of Gaussian variable, 
whence the proof follows. □ 
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In the case when 9 < 0 we can establish consistency, but not the strong 
consistency of (7) as n —>■ oo. According to Corollaries 3.3 and 4.2, we need to 
bound the following random variables: for 1 < m < jj 


Kn(m, a, /?) := 


1 mH-\-m—2H—l-\-ot 




r 

Jo 


X% ds + 


and K 2 := 




1 —1 y 

n 2^k =0 


k.n 


and for m > jj 


Kn(m, a, (3) := 


1 2Hm-\-m—2H—2-\-Oi 


r 

Jo 


X'i ds + dr 


and the same K 2 . 

Lemma 5.2. Let 9 < 0. 

(i) For any 1 < m < < -jj there exist such a > 0 and /3 > 0 that 

Kn(m, a, (3) ->• 0 

a.s. as n —» oo. 

(ii) For any Wry < m < jj there exist such a > 0 and /? > 0 that 


Kn(m, a, 0) 0 


in probability as n —>• oo. 


H 


(in) There exist such a > 0 and /3 > 0 that for any m > 

Kn(m,a,0) 0 


in probability as n —> oo. 

Proof, (i) Let 1 < m < < 77 - Then exponent mH+m — 2H — l is negative. 

Indeed, for Ff < | we have inequality 2 ^ f 1 1 < 77 . At first, choose a > 0 so that 
mFl + m — 2 FI — 1 + a < 0 and put $ = a. Then it is sufficient to note that 
n mH+m- 2 H-i+a _^ q an( j q a.s. as 7 i —»• oo while the integral f" X 2 ds 

is increasing with probability 1 and tends to nonzero with probability 1 random 
variable as n —> oo. 

(ii) Let y < in. < -7. Then exponent mFf + m — 2H — 1 is positive. Choose 
a = /3. It is sufficient to prove that there exists such a > 0 that 


—mH—m-\-2H-\-l—ot 


OO 


r-n 

: / X 2 s ds + n- mH - m+ 2 H+ 1 - a d n (a) 

J o 

in probability asn-> oo. In view of (34) it is equivalent to 

,n m ~ 1 

n -mH-m+2H+l-a / ds ^ ^ in probability as 77. — y OO. 
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To establish this convergence note that it follows from Cauchy-Schwarz in¬ 
equality that 





X 2 S ds > n -mH-2m+2H+2- a 



Denote 7 = —mH — 2m + 2 H + 2 — a < 0. Without loss of generality suppose 

ji rn ~ 1 

that *0 > 0. Note that f Q X s ds is a Gaussian process with mean 




X 0 ( Ori 

e \ 



and variance 


G 


2 

n 



EA'gA’t ds dt — e 2 n 


'f ‘ E ( ee "'I 


n —6u tdH 


Btf du + B 


H 


x 



„—6z tdH 


B“ dz + Ti¬ 


lt 


ds dt. 


(35) 


Since for any u,z > 0 we have E> 0, variance can be bounded from 
below by the following value 

/*n m—1 /*n m—1 

g^> / E B^B^dsdt 

Jo Jo 

= l n (m-l)(2H+2) f 1 r 1 ( S 2H + dt = Cn (m-l)(2H+2) _ 

2 Jo Jo 

(36) 

Note that other terms in (35) are of the same order so bound (36) is sharp. Now, 
denoting A/”(0,1) standard Gaussian random variable and 

$(:r) = ( 2 tt)“ 1 ? e _ T~ dy, we can deduce that for any A > 0 and sufficiently 
large n 


X s ds 1 < A 2 \ =~P In* 


= P 

< <& 


|nS \c7 n Af ( 0 , 1 ) + e„| < A| = $ ^ 


cr n n 2 &n 


$ _ 


cr n n2 
A e r . 


X s ds 


_ 


< A 
A 


< 


C 


(7 n ri2 CTr, 

C 


< 2 


tr„n2 

A 


< 


X — mH _ a 

G n nz n 2 2 



(37) 
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Choosing 0 < a < mH we get the proof of (ii). 

(Hi). For m > jj exponent 2Hm + m — 2 H — 2 is positive. Therefore, we 
repeat the proof of (ii) with the same a„ and with 7 = —2Hm — 2m+2H+3 — a 
instead of 7 . So, in the inequality similar to (36), we get in the right-hand side 
the upper bound 

C C 

cr n rff / 2 ~ 

Choosing 0 < a < \ we get the proof of (iii). 

□ 

Remark 2 . We can prove more than it was mentioned in (i), namely, to es¬ 
tablish that 

X 2 ds —> 00 
a.s. osn->oo (see Lemma 6.2 in Section 6). 

Lemma 5.3. Let m > ^jj. Then K 2 —>• 0 in probability as n —> 00. 

Proof. We apply the same method as in the proof of Lemma 5.2, but to the sum 
n ci " 1 X 2 n instead of integral J" Xf ds. As before, suppose that xq > 0. 
According to Cauchy-Schwarz inequality, 

.. n m —1 /n m —1 \ 2 

- E >«-”-■( E **,» , 

1 k—0 V fc=0 / 

where X]fc=o 1 fc,n is the Gaussian random variable with mean 



0 < e n = XQ ^2 - 


nx 0 


Os 1 / TIXq 

e ds < -— 


and variance 


k—0 


1 — ^2 ^Xk,nXj^n > ^2 
k,j =0 k,j =0 

A n m — 1 

I V ((±) 2H + (±) 2H - |i - t\ 2H ) 

9 / j y\nJ ' \nJ In nl J 

k,j =0 

n m — 1 n m — 

E E 




= - Z. (fJ ^ m ~i) 

fc=0 Z=0 

n m —1 n m —1 

-2H+m(2H+l)+m t l \ 2H+1 1 


= n" 2ff £ Z 2ff+1 = 


n 


E (^r 


2=0 


z=o 


^-j n —2H+m(2H+l)+m 
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Therefore, for any e > 0 and x n > 0 


ICnl n m 2H _ I , „ I 1 


n — 1 


P{Kl>e)<P -JSlL__> e < p - £ Xt, n < 


E n ""-1 ^2 
/c—0 


/c—0 


+P{ICn| >^en 2if “ m } <P 


1 <„L. E «- 


1 m+l 


k—0 


< p 


Similarly to (37), 


nm+l 


^ ^ Xk,n I ^ %n / + 


c 


fc=o 


1 


^ra+l 


E X ^ 


< X n > = P 


k—0 


jlo-n^O, 1) + e n \ < xZn 7 "* 1 1 


- $ [ _ 2 hl + ] _ $ (_£hl _ x ^ nTn ^ 1 

\ O n O n I \ O n O n 


< c 


' ^ 22+1 
Xn n 2 n 


Evidently, for any m > 1, — —> 0 as n — > oo. Therefore, to supply con- 

On 

vergence K 2 —>• 0 in probability, we need to choose x n in such a way that 

I m +1 

1 rill 2 

—— , „ „ — —» 0 and ——- —>• 0 as n —> oo. Put x n = n r . Then r must 

/y>Z rv-i 4il Zm sy 

•^n 1 u n 

satisfy the double inequality 

m — 2 H < r < m + 2 Hm — 2 H — 1. 

This inequality can be satisfied only for m > whence the proof follows. □ 

Theorem 5.4. Let m > Then the estimator 9 n (m), introduced in (7), is 
consistent. 


6. Auxiliary results 

At first we establish an auxiliary result concerning the bounds for several sums 
of integral type that will participate in the bounds for the numerator of (7) 

Lemma 6.1. For any m > 1 and n> 2 there exists C > 0 not depending on n 
such that 

(i) 

n m _l 

E H log 2 < Cn^ m ~ 1)H+rn log 2 n, 

k=0 
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(a) 

n m — 1 

53 (*±i ) 2H log 4 log 4 n. 

k—0 

Proof. We base the proof of both statements on the following evident relation: 
for any function / : [0,1] —>• R that is Riemann integrable on [0,1], and for 
any m > 1 integral sums S(f(x),n m ) := -^rJ2k=o 1 f tenc ^ to integral 

fo f{ x )dx as n —> oo. In particular, these integral sums are bounded. Consider 
the statement (i). Evidently, 

n m — 1 n m —1 

£ (i±.)"log^ =„<-«« •£ (^"log^n- 1 ) 

k—0 k—0 

n m — 1 

<2 n^ H ^ (^)"(log 2 (^)+(m-l) 2 log 2 n) 
k—0 

= 2 n {m - l)H+m S(x H log 2 x, n m ) + 2(m - 1) log 2 n • 5 (x 44 , n m ) 

< (7n ( ”- 1)fl+ra + Cn( m - 1)4f+m log 2 n < Cn (m " 1)44+m log 2 n 


for n > 2. Statement (ii) is established similarly. 

□ 

Next auxiliary result establishes asymptotic behavior of integral jf Xf ds as 
T —f oo. 

Lemma 6.2. Let process X satisfy equation (1). Then ff Xf ds —»• oo with 
probability 1 as T —>• oo. 


Proof. The result is obvious for 9 > 0 therefore we consider only the case 
9 < 0. Since jf X^ ds is nondecreasing in T > 0, it is sufficient to prove that 
ff Xf ds —» oo in probability. For any A > 0 consider the moment generation 
function 0p(A) = Eexp|— A ff X% dsj and 0oo(A) = Eexp{— A / 0 °°X 2 ds} 
so that 

0oo(A) = lim 0 r (A). 

T—too 

Evidently, 


whence 


X 2 ds > T" 1 



0 t(A) < 0 ^(A) 



Random variable T 2 ff X s ds is Gaussian with mean m(T ) and variance 
ct 2 (T), say. Note that for a Gaussian random variable £ = m + <rAf( 0,1) we 
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have that 

Eexp {-AC 2 } = (2Aa 2 + l)~' 2 exp |- + 1 | < (2A<r 2 + l)~U 

Therefore, it is sufficient to prove that 

lim <j 2 (T) = oo. 

T-> oo 

Similarly to (36). 

o*(T) > T 2H+1 £ £ (s 2H + t 2H -\s- t\ 2H ) dsdt^oo 
as T —> oo, whence the proof follows. □ 


7. Simulations 

In this section, we present the results of simulation experiments. We simulate 
20 trajectories of the fractional Ornstein-Uhlenbeck process (1) with .To = 1 
for different values of 9 and H. Then we compute values of 9 n {m). For each 
combination of 0, H, n and m the mean of the estimator is reported. 

In Tables 1-3 the true value of the drift parameter 9 equals 2. In this case 
the behavior of the estimators is almost the same for different values of H. Also 
we can see that the value of 9 n (m) is determined by n and does not depend on 
m. Further, we consider the case of negative 9. We simulate the process with 
H = 0.45, 9 = —3 and m = 4, 5. The results are reported in Tables 4-5. One 
can see that the method works but the rate of convergence to the true value 
of a parameter is not very high. There are two reasons for this: the estimator 
is only consistent not strongly consistent and moreover, the trajectories are so 
irregular that even the length of the interval is small we can not “catch” the 
trajectory. 


Table 1 
e = 2, m = 2. 


n 

5 

10 

50 

100 

500 

1000 

H = 0.05 

2.45763 

2.21281 

2.0395 

2.01911 

2.00300 

2.00100 

H = 0.25 

2.45766 

2.21281 

2.0395 

2.01911 

2.00300 

2.00100 

H = 0.45 

2.45794 

2.21281 

2.0395 

2.01911 

2.00300 

2.00100 
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